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La presente invention conceme une methode pour obtenir des temps de trajet en 
reflexion a partir d'une interpretation de donnees sismiques en ondes cylindriques migrees, 
pour une valeur donnee d'un parametre definissant la pente de ces ondes, ou de la 
5 superposition de telles donnees associees a differentes lignes d' acquisition sensiblement 
paralleles, ce parametre pouvant eventuellement prendre successivement plusieurs valeurs. 

ETAT DE LA TECHNIQUE 

La demigration de donnees sismiques est I'operation inverse de la migration qui, on 
le rappelle, est une methode classique d'imagerie de donnees sismiques, consistant 
essentiellement, connaissant la valeur d'un champ d' ondes a une profondeur connue, par 
exemple en surface, ainsi qu'un modele de la repartition des vitesses de propagation des 
ondes dans le sous-sol, a modeliser la propagation du champ source et la retropropagation 
des donnees de reflexion enregistrees et a chercher des coherences de phase entre ces deux 
champs modelises. La migration est particulierement utile pour interpreter des donnees 
sismiques acquises sur des structures complexes. 

EUe permet aussi d'acceder, via une demigration, aux temps d'arrivee des reflexions 
associees aux evenements pointes. Cette operation intervient, notamment, lorsqu'il s'agit de 
determiner, par des methodes cinematiques telles que les methodes d'analyse de vitesse de 
migration, la distribution de vitesse dans la subsurface. Cette determination constitue une 
20 etape cle dans I'imagerie des structures complexes. Un exemple d'une telle methode est 
decrit par : 

Stork, C, 1992 ; « Reflection Tomography in the Postmigrated Domain ; Geophysics, 
57,680-692. 

Une methode d'analyse de vitesse de migration particulierement performante dite 
25 Smart (pour « Sequential Migration Aided Reflection Tomography ») est decrite 
notamment par : 

- Ehinger, A. et Lailly P., 1995, Velocity Model Determination by the SMART 
Method, Part 1 : Theory : 65^*^ Ann. Intemat. Mtg., Soc. Expl. Geophys., Expanded 
Abstracts, pages 739-742. 
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C'est a la tomographic de reflexion qu'est confiee, dans les methodes de ce type, la 
tache de determiner le modele de vitesse ; les temps de trajet en reflexion (donnees de la 
tomographic de reflexion) sont obtenus, dans les parties complexes, en demigrant les 
evenements pointes lors de I'interpretation des donnees migrees. Cette methode a ete mise 
en oeuvre avec succes pour Timagerie en 3D de structures complexes en utilisant des 
algorithmes de migration 2D avant sommation (et la demigration associee) sur un 
ensemble de prof i Is sismiques extraits d'une acquisition 3D. Neanmoins une telle approche 
ne s'avere efficace que dans un contexte bien specifique : 

sur le plan geologique la structure doit varier relativement lentement suivant une 
direction d*espace ( appelee « strike » dans la litterature anglo-saxonne) ; 

- les profils sismiques doivent etre acquis suivant la direction orthogonale (done la 
direction de complexite) appelee "Mip" dans la litterature anglo-saxonne. 

En revanche, cette approche peut difficilement etre envisagee sur des structures 
complexes sans direction privilegiee : sur de telles structures il faut avoir recours a des 
techniques de migration 3D avant sommation, techniques gourmandes en temps de calcul 
surtout si rimagerie doit etre realisee sur un volume. 

Par le brevet FR 2 784 195 du demandeur notamment, on connait une autre 

methode de migration de donnees sismiques, connue des gens de Fart sous le nom 
d'Alchemig. Cette methode permet de realiser en 3D, une migration avant sommation 
d'evenements sismiques, et d' imager des volumes d'une zone souterraine, ceci pour un 
temps de calcul attrayant. a partir d'une serie d'un nombre Ns de cycles de sismique 
reflexion. 

Elle comprend remission successive de champs d'onde elementaires defini chacun 
par Tassociation d'un signal sismique W(t) et d'un lieu d'emission defini dans une serie de 
— > 

lieux d'emission Si avec 1 < i < Ns, la reception par des recepteurs sismiques places en des 

positions R ^ , des signaux sismiques renvoyes par la zone en reponse a chacun de ces 
i 

champs d'onde, et I'enregistrement des differents signaux regus par chaque recepteur 



sismique sous la forme de traces sismiques d ^ (t) dependant du temps. Pour un modele de 

i 

vitesses donne. 



a) on definit un vecteur de lenteur p (homogene a T inverse d'une vitesse) dont les 
deux composantes py et px peuvent chacune prendre une suite de valeurs prealablement 

5 definie, 

— > — » 

b) on definit pour un vecteur de lenteur p donne et lieu d' emission Si donne, une 

fonction de decalage temporel to( p , i) ; 

c) on applique la fonction de decalage temporel to( p , i) a chaque champ d'onde 

elementaire associe au lieu d'emission Siet on forme un champ d'ondes composite en 
10 surface par superposition spatio-temporelle des differents champs d'onde elementaires 
ainsi decales ; 

j , , 

d) on applique un decalage temporel to( p , i) a chaque trace sismique d (t) reperee 

i 

par le couple (i,j) et on forme un champ de traces composite en surface par superposition 
spatio-temporelle des differentes traces sismiques ainsi decalees ; 

15 e)on effectue une migration du champ de traces composite en utilisant comme 

champ d'ondes le champ d'ondes composite, ceci en modelisant la propagation de champ 
d'ondes composite ainsi que la retropropagation de champ de traces composite, et en 
combinant de fa5on adaptee les deux champs composites ainsi modelises en tout point de 
la zone a imager ; 

20 f) on repete les etapes c) a e) pour toutes les valeurs prises par les composantes py et 

Px du vecteur p ; et 



g)pour toute valeur fixee de la deuxieme composante px du dit vecteur p , on 
somme les resultats de ces differentes combinaisons de maniere a obtenir une image 
migree associee a cette valeur fixee de px, realisant ainsi une migration avant sommation. 




On suppose, meme si la technique s'avere robuste en cas de violation de ces 
hypotheses, que les donnees sismiques sont a azimut constant, (razimut etant defini 
comme la direction des bipoints source-capteur) (hypothese 1), et que I'acquisition est 
effectuee en deplagant la source suivant des lignes paralleles a la direction definie par 
5 I'azimut (hypothese 2). 

Si ce n'est pas le cas, on peut s'y ramener en procedant dans un premier temps a un 
pretraitement suivant par exemple la technique dite AMO pour « Azimut Moveout, » bien 
connues des gens de I'art, pour amenager les donnees provenant d'acquisitions marines 
standards. 

10 Le systeme de coordonnees choisi ici, est tel que Taxe des x est parallele a la 

direction des lignes d' acquisition et Taxe des y repere done la position d'une ligne 
d' acquisition. Ci-apres, nous definirons par « donnees en ondes cylindriques » la mesure, 
suivant une ligne d'acquisition, de la reponse sismique de la subsurface a une onde 
cylindrique (onde engendree par une ligne de sources et dont la phase varie lineairement 

15 avec Tabscisse suivant la ligne d'acquisition) dont I'axe coincide avec la ligne d'acquisition 
consideree. Compte tenu de cette definition, la technique Alchemig permet d'obtenir la 
superposition, suivant les differents lignes d'acquisition, de « donnees en ondes 
cylindriques migrees », ces ondes cylindriques etant associees a une ligne d'acquisition et a 
une valeur predefinie du parametre px definissant la pente de Tonde cylindrique. La ^i 

20 demarche consiste a appliquer aux donnees des dephasages lineaires parametres par le ^ 

vecteur ^dont la composante suivant I'axe des profils est px, et a calculer une telle 
superposition en appliquant aux donnees dephasees un logiciel de migration par ondes 
planes . La sommation apres migration s'obtient en sonmiant les resultats obtenus pour les 
differentes valeurs prises par px. On suppose, par un choix convenable de Torigine du 
25 repere ou de la fonction de dephasage, que le point de la surface ou les dephasages sont 
nuls a pour coordonnees (x = 0,y = 0). 

On note que la technique presentee ci-dessus, constitue une fagon parmi d'autres 
d'obtenir la superposition de donnees en ondes cylindriques migrees associees a des lignes 
d'acquisition paralleles et a une valeur donnee du parametre px- Une autre fagon de 
30 proceder peut consister par exemple, comme decrit par Claerbout, J.F.,1971, « Tov^ards a 
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Unified Theory of Reflector Mapping ; in Geophysics, 36, N°3, 467-481, a calculer le 
resultat en effectuant pour chaque ligne d'acquisition : 

la propagation de Tonde cylindrique et la retropropagation des donnees en ondes 
cylindriques ; et 

5 - le calcul, en tout point de Tespace, de rtntercorrelation du champ source dont on a 
modelise la propagation et du champ de donnees retropropagees, ce qui permet 
d'obtenir, pour la ligne d*acquisition consideree, les donnees en ondes cylindriques 

migrees ; 

et en sommant (superposant) les donnees en ondes cylindriques migrees obtenues pour les 
10 differentes lignes d'acquisition. 

Si Ton dispose, pour plusieurs valeurs du parametre px, de la superposition de 
donnees en ondes cylindriques migrees associees a differentes lignes d'acquisition, on peut 
comme explique dans le brevet precite, controler la qualite du modele de vitesse utilise 
pour realiser la migration et aussi mettre a jour le modele de vitesse en exploitant les 
15 deformations des evenements lorsque Ton passe d'un parametre d'onde cylindrique a un 
autre (analyse de vitesse de migration). 

Pour proceder a cette mise a jour, comme il est fait suivant la methode dite Smart 
precitee (ou d'autres), il faut, pour differentes valeurs du parametre px. pointer les 
evenements dans la superposition des donnees en ondes cylindriques migrees, 
20 superposition correspondant aux differentes lignes d'acquisition, puis demigrer les 
evenements pointes et confier a la tomographic de reflexion le soin de la mise a jour du 
modele. Encore faut-il savoir comment demigrer les evenements migres interpretes lorsque 
r interpretation a ete effectuee sur des donnees en ondes cylindriques migrees ou leur 
superposition obtenue pour differentes lignes d'acquisition. 

25 La methode selon I'invention 

C'est un des problemes complexes que la methode selon I'invention, se propose de 
resoudre. L'approche retenue debouche, de plus, sur des methodes que Ton decrira par 
ailleurs, qui sont connexes au concept d'onde cylindrique. 




La methode selon T invention permet de determiner les temps de trajet en reflexion 
d'evenements sismiques pointes sur des enregistrements en 3D de donnees sismiques 
correspondant a des signaux captes par des recepteurs sismiques repartis suivant une ligne 
d'acquisition (respectivement suivant plusieurs lignes d' acquisition) en reponse a 
remission dans le sous-sol d'ondes depuis un ou plusieurs points-source, ces donnees 
ayant ete au prealable transformees en donnees en ondes cylindriques migrees et 
interpretees (respectivement en une superposition donnees en ondes cylindriques migrees 
et interpretees). 

Dans le cas oii 1' interpretation a ete effectuee sur les donnees migrees associees a 
une seule onde cylindrique, et done une seule ligne d' acquisition, la methode comporte les 
etapes suivantes : 

a) on definit un vecteur lenteur ( p ) dont la composante (px) suivant une 
direction parallele a la ligne d' acquisition, definit la pente de I'onde cylindrique ; 

b) pour un recepteur sismique situe a I'abscisse ( a:^ ) sur la ligne d'acquisition, 
on recherche I'abscisse (<^) du point-source sur la ligne d'acquisition telle qu'un rayon 
partant du recepteur sismique et se reflechissant sur un evenement pointe, emerge au dit 

point source, avec un vecteur lenteur ( p ) dont la composante suivant le ligne d'acquisition 
est (-px) ; 

c) on determine un temps de trajet ( (Xj^ ) ) en ajoutant a la valeur du temps de 
trajet le long du rayon, un temps egal au produit de la pente de I'onde cylindrique par 
Tabscisse du point source ; 

d) on repete les etapes b) et c) pour toutes les positions des recepteurs pour 
lesquelles on veut un resultat de demigration ; et 

e) on repete les etapes a) a d) pour toutes les valeurs prises par le parametre 
(Px) pour lesquelles on a pointe un evenement dans les donnees migrees correspondantes. 

Dans le cas ou T interpretation a ete effectuee sur la superposition de donnees en 
ondes cylindriques, superposition realisee a partir de plusieurs lignes d'acquisition, la 
methode comporte les etapes suivantes : 



a) on definit un vecteur lenteur ( p ) dont la composante (px) suivant une 
direction parallele aux lignes d' acquisition, definit la pente des ondes cylindriques 
associees aux differentes lignes d' acquisition ; 

b) pour un recepteur sismique donnee situe a I'abscisse (x^) sur une ligne 
d' acquisition, on recherche I'abscisse ) d'un point source sur la ligne d'acquisition telle 

qu'un rayon partant du recepteur sismique et se reflechissant sur un evenement pointe, 

— > 

emerge au dit point source, avec un vecteur lenteur ( p ) dont la composante suivant la 
ligne d'acquisition est (-px) ; 

c) on determine un temps de trajet (r' (;c^ ) ) en ajoutant a la valeur du temps de 
trajet le long du dit rayon, un temps egal au produit de la pente de I'onde cylindrique par 
I'abscisse du point source ; 

d) on repete les etapes b) et c) pour toutes les positions des recepteurs pour 
lesquelles on veut un resultat de demigration ; et 

e) on repete les etapes a) a d) pour toutes les lignes d'acquisition pour 
lesquelles on veut un resultat de demigration ; 

f) on repete les etapes a) a e) pour toutes les valeurs prises par le parametre 
(Px) pour lesquelles on a pointe un evenement dans les donnees migrees correspondantes. 

Le traitement prealable aboutissant a des donnees en ondes cylindriques migrees, 
peut eventuellement inclure si necessaire un pretraitement de type AMO deja cite ou un 
reechantillonnage spatiales des donnees sismiques acquises sur le terrain. 

Suivant un mode de mise en ceuvre, on utilise les temps de trajet associes aux 
evenements en ondes cylindriques, pour mettre en oeuvre une technique de tomographic de 
reflexion par ondes cylindriques et calculer une distribution des vitesses dans le milieu. 

Suivant un autre mode de mise en ceuvre, on convertit les temps de trajet associes 
aux evenements en ondes cylindriques, en temps de trajet associes aux evenements par 
points de tir, en exploitant les correspondances entre les donnees en ondes cylindriques et 
les donnees en points de tir. 
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Suivant un mode de mise en oeuvre, lorsque le dispositif d'acquisition des donnees 
sismiques est organise par lignes d' acquisition, c'est-a-dire lorsque sur une meme ligne 
figurent tout un ensemble de points-source et tout un ensemble de capteurs, on transforme 
les temps de trajet par tir synthetise en des temps de trajet associes aux couples source - 
recepteurs (du dispositif d' acquisition des donnees sismiques), cette transformation se 
faisant par extrapolation suivant la direction des lignes d' acquisition, du temps de trajet 
synthetise pour le couple source - recepteur le plus proche du couple source - recepteur 
que Ton veut documenter en temps de trajet. 

Suivant un mode de mise en oeuvre, les donnees sismiques etant seuelement a 
azimut constant, on transforme les temps de trajet par tir synthetises, en des temps de trajet 
associes aux couples source - recepteurs cette transformation se faisant par une double 
extrapolation, la premiere suivant la direction definie par I'azimut, la deuxieme suivant la 
direction orthogonale a la direction precedente. 

La methode peut etre utilisee notanrmieht pour la determination du modele de 
Vitesse en appliquant par exemple aux temps de trajet obtenus, une tomographie de 
reflexion. 

PRESENTATION SOMMAIRE DES FIGURES 

D'autres caracteristiques et avantages de la methode selon I'invention, apparaitront a 
la lecture de la description ci-apres d'exemples non limitatifs de mise en ceuvre, en se 
referant aux dessins annexes oii : 

la Fig.l donne un descriptif des differentes options permettant la mise en oeuvre de la 
methode selon T invention ; et 

la Fig.2 confronte, sur une section sismique a offset constant, les donnees sismiques 
avec les temps d'arrivee obtenus par mise en oeuvre d'une option de la methode. 

DESCRIPTION DETAILLEE 

La demigration d'evenements pointes sur des donnees en ondes cylindriques 
migrees, demande une technique specifique. Nous distinguerons successivement : 

la demigration d'un evenement pointe sur des donnees en onjdes cylindriques migrees ; 




la defhigration de la superposition (« stack ») de donnees en ondes cylindriques 
migrees, associees a differentes lignes d'acquisition ; 

Tobtention de donnees cinematiques (temps de trajet en reflexion)associees a des 
couples source-capteur ; et 

5 - Tobtention de donnees cinematiques associees a des couples source-capteurs existant 
dans le dispositif d*acquisition d'origine (c'est~a-dire les couples source-capteurs 
utilises pour realiser la migration). 

1) Demigration d'un evenement pointe sur des donnees en ondes cvlindriques 

migrees 

10 Le resultat d'une telle demigration est le temps d'arrivee, sur les enregistrements 

associes a Tonde cylindrique consideree (caracterisde par la coordonnee suivant Taxe des y 
de la ligne d'acquisition et par la valeur du parametre px definissant la pente de Tonde 
cylindrique), de la reflexion correspondant a I'evenement pointe. Ce temps d'arrivee r*' 
depend de I'abscisse .r^ du recepteur considere. La demigration s'effectue naturellement 

15 en recherchant, pour une position de recepteur donnee, 1' abscisse ^ du point source ou des 

points sources 5 sur la ligne consideree, tel que le rayon partant du recepteur et se 

— > 

reflechissant sur I'evenement pointe emerge en S avec un vecteur lenteur p dont la 
composante suivant x vaut -px, le modele de vitesse etant identique a celui utilise pour 
realiser la migration. 

20 Ainsi formule, le probleme apparait comme la resolution d'une equation non 

lineaire (exprimant la condition sur le vecteur lenteur) a une inconnue (I'abscisse £ de la 
source). Cette equation peut etre resolue numeriquement par n'importe quelle methode 
adaptee a la resolution d'equations non lineaires telle que la methode de Newton par 
exemple. 

25 Le temps de trajet resultat de la demigration t^'iXj^) sera obtenu en prenant la 

valeur du temps de trajet le long de tous les rayons obtenu pour la source S obtenu pour la 
source S, temps auquel on aura ajoute la quantite p^^ . On note qu41 y a, en fait, une 

infinite de tels rayons mais ces rayons foumissent, pour une arrivee donnee, le meme 
temps de trajet. L'operation est repetee pour toutes les abscisses de recepteurs pour 



lesquelles on souhaite avoir le resultat de demigration. Les temps d'arrivee ainsi obtenus 
constituent une information intrinsequement liee aux donnees en ondes cylindriques 
synthetisees. lis ne dependent pas du modele de vitesse utilise pour effectuer la rriigration, 
(ce qui justifie qu'on puisse les utiliser pour la determination du modele de vitesse) car 
reoperation de demigration presentee ci-dessus est, sur le plan de la cinematique, I'inverse 
de la migration de donnees en ondes cylindriques. La composante suivant x du vecteur 

lenteur p du rayon a I'emergence au capteur sismique constitue egalement une quantite 
intrinsequement liee aux donnees et qui ne depend done pas du modele de vitesse utilise 
pour la mise en oeuvre de la migration ; elle donne la pente, sur les donnees en ondes 
cylindriques, de Tevenement enregistre au capteur, 

2) Demigration de la superposition de donnees en ondes cylindriques migrees. 
associees a differentes lignes d 'acquisition 

Cette operation s'effectue simplement en repetant I'operation definie au paragraphe 
precedent pour toutes les lignes d'acquisition ou Ton souhaite obtenir un resultat de 
demigration et elle delivre les memes quantites intrinseques, a savoir le temps d'arrivee de 
Tevenement, sur les donnees en ondes cylindriques associees a chaque ligne, en fonction de 
Tabscisse du capteur ainsi que la pente de cet evenement suivant I'axe des x, c'est a dire la 
composante suivant x du vecteur lenteur lors de I'emergence du rayon au capteur. De 
plus, si nous appelons pi et pi les composantes suivant I'axe des y du vecteur lenteur 
pour des rayons emergeant en surface respectivement au capteur et a la source, la quantite 
Pr + Ps apparait elle aussi comme une quantite intrinseque ; elle s'interprete comme la 
pente, suivant I'axe des y (et done quand on passe d'une ligne a I'autre), de Tevenement vu 
pour I'abscisse de capteur sismique considere, sur des donnees en ondes cylindriques 
associees au parametre px. 

3) Obtention de donnees cinematiques (temps de traiet en reflexion) associees a 
des couples source-capteurs 

Si Ton dispose d'un logiciel de tomographie de reflexion pouvant traiter des temps 
de trajet associes a des ondes cylindriques, les donnees obtenues au paragraphe precedent 
suffisent pour permettre la mise a jour du modele de vitesse et realiser ainsi, suivant la 
methode dite Smart, I'analyse de vitesse de migration. 



Si Ton ne dispose pas de tels logiciels d'inversion, mais seulement de logiciels 
classiques manipulant des donnees cinematiques organisees par couple source-recepteur, 
on realisera I'etape decrite ci-apres. 

Cette etape consiste a exploiter des correspondances entre les donnees en ondes 
cylindriques et les donnees par point de tir, c'est-a-dire Tensemble des enregistrements 
effectues par les differents recepteurs pour une menne source sismique supposee 
ponctuelle. 

Une telle correspondance ne peut etre exploitee que si Ton dispose de donnees en 
ondes cylindriques associees a plusieurs valeurs du parametre definissant la pente des 
ondes cylindriques. Notons que la migration ne joue ici aucun role : il s'agit simplement de 
transformer des donnees. Les donnees cinematiques (temps de txajet en reflexion) associees 
a des couples source-capteur sont obtenues en procedant ligne d'acquisition par ligne 
d'acquisition (les couples source-capteur sont, rappelons-le, supposes situes sur ces lignes). 

Pour une ligne donnee, on dispose, par application des techniques presentees au 
deux paragraphes precedents, de la fonction r''(jc^), ceci pour differentes valeurs du 
parametre . En d'autres termes, nous disposons de la fonction, que nous designons, au 
prix d'un changement de notation, t'^ip^.x^) ou, plus exactement, d'un echantillonnage de 
cette fonction, les valeurs prises par cette fonction ayant ete obtenues pour des valeurs 
specifiques du couple {p ^ ,x^). On definit une fonction ip^^.^R ) par la relation : 

^S^Px^X^) = -—{p„Xi,) (1) 

quantite que nous pouvons evaluer numeriquement (par exemple en utilisant des 
approximations basees sur des differences finies) pour des valeurs specifiques du couple 
(,p^,x^). La formule (2) ci-apres nous donne le temps note t"{x^,Xn) de I'dvenement au 
capteur d'abscisse x^ dans les donnees par point de tir pour une abscisse de tir 



T'{Xs,x^) = t'{p^,x^)- p^Xs(,Px,x^) 



(2) 




Les temps ainsi obtenus sont des quantites intrinseques : ils sont uniquement lies 
aux donnees sismiques et ne dependent pas de la fagon suivant laquelle la migration a ete 
mise en oeuvre. En conclusion, il suffit d'appliquer la formula (2) sur les differentes lignes 
pour tous les couples (/?^,a-^) pour lesquels on dispose des donnees t'^iPj^.-^n) et 

5 (p ,Xf^) pour connaitre les temps d'arrivee de l*evenement interprete (ou, dans la 

pratique, de chaque evenement interprete) pour une multitude de couples source-capteur 
couvrant la surface et permettre ainsi une detennination precise, via une tomographie de 
temps de trajet en reflexion, du modele de vitesse. 

4) Obtention de donnees cinematiques associees a des couples source-capteur 
10 existant dans le dispositif d'acquisition d'origine 

Meme si les donnees obtenues par la technique exposee a la section precedente 
suffisent pour effectuer la mise a jour du modele de vitesse, il peut s'averer utile de 
confronter les temps obtenus avec les donnees sismiques d'origine non migrees. Une telle 
confrontation est en effet interessante pour valider le pointe (et notamment Tinterpretation) 
15 effectue sur les donnees migrees et pour s'assurer de la fiabilite de la demigration qui, on 
I'a vu, fait appel a des procedures numeriques basees sur des approximations ainsi qu'au 
trace de rayons, outil tres sensible en la representation numerique des reflecteurs. 

Cependant une confrontation directe peut difficilement etre envisagee car les 
couples source-capteur n'ont aucune raison d'etre ceux correspondant aux donnees utilisees 

20 pour realiser la migration. En effet, les temps de trajet foumis par la methode exposee a la 
section precedente, sont obtenus pour des couples source-capteur repartis suivant une suite 
de nuages de points distribues dans des espaces 2D (chacun de ces espaces correspondant a 
une ligne d'acquisition) alors que les donnees d'origine peuvent etre repartis suivant un 
nuage de points repartis dans un espace 3D. Cependant meme dans le cas ou la migration a 

25 ete realisee a partir de donnees organisees par lignes 3D, la difficulte demeure : le couple 
source-capteur defini par le couple d'abscisses ( X^(p^,Xf^),Xf^) , obtenu comme expose a 
la section precedente, n*a aucune raison de figurer dans le dispositif d'acquisition d'origine. 

Pour obtenir des donnees cinematiques associees a des couples" source-capteur 
existant dans le dispositif d'acquisition d'origine a partir des resultats obtenus suivant la 
30 methode exposee a la section precedente, il faut avoir recours a des techniques 



d'interpohition ou d'extrapolation. Une methode consiste a se definir un critere de 
proximite dans un espace 2D ou 3D suivant le cas et, pour un couple source-recepteur 
donne du dispositif d'acquisition d'origine, y affecter le temps de trajet du couple source- 
capteur le plus proche parmi ceux qui sent documente en temps de trajet par mise en oeuvre 
de la methode exposee a la section 3. Une telle approche est viable si la distance 
d^extrapolation de Tinformation n'est pas trop grande. Dans le cas contraire il faut avoir 
recours a des techniques plus elaborees telles que celles exposees ci apres, selon que Ton 
opere sur des donnees reechantillonees ou sur des donnes a azimut constant. 

4-1 Mise en ceuvre sur des donnees organisees par lignes d 'acquisition 

Dans ce type de mise en oeuvre, nous supposons aussi que les donnees sismiques 
sont a azimut constant et que I'acquisition est effectuee en deplagant la source suivant des 
lignes paralleles a la direction definie par Tazimut suite par exemple a un reechantillonnage 
spatial des donnees et/ou un pretraitement de type AMO, conmie on I'a vu. Les couples 
source-capteur intervenant dans la mise en oeuvre de la methode de migration sont repartis 
suivant une suite de nuages de points distribues dans des espaces 2D, chacun de ces 
espaces correspondant a une ligne d'acquisition. 

Dans ce cas, les donnees sont organisees par ligne d'acquisition et nous voulons 
obtenir le temps d'arrivee de I'evenement reflechi pour le couple source-capteur du 
dispositif d'acquisition d'origine (disp), defini par la ligne d'acquisition a laquelle il 
appartient et par les abscisses jc/"'' et x^"^ . Nous avons, comme explique ci-dessus, extrait 
de Tensemble des temps d'arrivee calcules pour la ligne consideree, celui associe au couple 
source-capteur defini par le couple (a*^"''' , jc^f'' ) qui s'averait le plus proche de celui associe 

au couple ( xf'^ , x^"'' ) parmi le nuage de points obtenus a I'interieur de Tespace 2D decrit 
par les couples (x^ ,.y^ ) . Nous utilisons alors la formule (3) qui suit, pour obtenir le temps 
d'arrivee t*' {xf'^ ,x^"^) a partir des quantites rendues disponibles au cours des eiapes de 
traitement precedentes : 

(xf' , ) = t' {xf^ , xf^ ) - {x'r - ^xf' ) + Pr (4"' - ) (3) 

la quantite ayant ete definie aux sections 1, 2 ci-dessus. On remarque que cette 
formule exploite le resultat important suivant : sous les hypotheses 1 et 2 enoncees a la 




section I, pour un capteur donne et done pour une ligne d'acquisition donnee, la pente d'un 
evenement dans les donnees en ondes cylindriques associees au parametre p^, est 
identique a la pente de I'evenement dans la collection point de tir pour le tir lie par la 
relation (1) au couple (p,,a-/j) a Tabscisse (.v^) du capteur. La technique d'extrapolation 
indiquee ci-dessus, n'est bien entendu pas limitative. On peut tout aussi bien utiliser une 
technique d' interpolation. 

4-2 Mise en oeuvre de la methode sur des donnees a azimut constant 

Dans ce type de mise en ceuvre, nous supposons que les donnees sismiques sont a 
azimut constant. Les couples source-capteur intervenant lors de la mise en ceuvre de la 
migration, sont repartis suivant un nuage de points distribues dans un espace 3D, ces 
couples pouvant 6tre decrits par Tintermediaire de trois parametres : les abscisses 
respectives de la source et du capteur et leur ordonnee commune. Le probleme consiste ici 
a obtenir les temps d'arrivee de Tevenenient reflechi pour les couples source-capteur du 
dispositif d'acquisition d'origine (distribues done suivant un nuage de points dans un 
espace 3D) a partir des temps calcules en utilisant la methode presentee a la section 3. 
Cette methode est alors mise en oeuvre pour une suite de lignes d'acquisition fictives 
reperees par leuf ordonnee y.,i, i representant le numero, de la ligne consideree. Ces 
lignes doivent etre suffisamment rapprochees pour permettre Textrapolation, sur les 
couples source-capteur du dispositif d'acquisition d'origine, des temps de trajet calcules 
suivant la methode presentee ci-dessus a la section 3. Nous voulons done calculer le temps 
de trajet de I'evenement reflechi pour un couple source-capteur du dispositif d'origine, ce 
couple etant repere par le couple d'abseisses ( xf'^ , x^^^^ ) et par leur ordonnee commune 

ydiSp 

Pour ce faire, nous extrayons de Tensemble des temps d'arrivee calcules pour la 
suite de lignes mentionnee ci-dessus, eelui associe au couple source-capteur defmi par le 
couple Cr^''''' , -v^'''' ) et I'ordonnee y. qui s'avere le plus proche, au sens de la metrique 
definie dans Tespace 3D, de celui associe au couple source-capteur considere du dispositif 
d'origine. Nous utilisons alors la formule (4) ci-apres pour obtenir le temps d'arrivee 
t^'ixf^x'^^^.y'^''^) a partir des quantites rendues disponibles au cours des etapes de 
traitement preeedentes. 
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La encore, la technique d'extrapolation est indiquee a titre d'exemple non limitatif. 
On peut tout aussi bien utiliser a la place, une technique d' interpolation. 

La figure 2 confronte, sur une section a deport (offset) constant (en Toccurrence 
1500 m), obtenue a la suite d' operations de prospection sismique 3D en Mer du Nord, les 
donnees sismiques avec les temps d'arrivee obtenus par mise en oeuvre de la technique 
exposee ci-dessus. On y retrouve clairement la reflexion sur la base de la structure salifere. 
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REVENDICATIONS 

1) Methode pour determiner les temps de trajet en reflexion d'evenements 
sismiques pointes sur des enregistrements en 3D de donnees sismiques correspondant a des 

5 signaux captes par des recepteurs sismiques repartis suivant une ligne d' acquisition en 
reponse a remission dans le sous-sol d'ondes depuis des points-source, ces donnees ayant 
ete au prealable transformees en donnees en ondes cylindriques migrees et interpretees, 
caracterisee en ce que : 

a) on definit un vecteur lenteur ( p ) dont la composante (px) suivant une 
10 direction parallele a la ligne d' acquisition, definit la pente de I'onde cylindrique ; 

b) pour un recepteur sismique situe a I'abscisse (.x^) sur la ligne d' acquisition, 
on recherche I'abscisse (^) du point source sur la ligne d'acquisition telle qu'un rayon 
partant du recepteur sismique et se reflechissant sur un evenement pointe, emerge au point 

source, avec un vecteur lenteur ( p ) dont la composante suivant le ligne d'acquisition est (- 
15 Px) ; 

c) on determine un temps de trajet ( r (jc^ ) ) en ajoutant a la valeur du temps de 
trajet le long du dit rayon, un temps egal au produit de la pente de I'onde cylindrique par 
I'abscisse du point source ; 

d) on repete les etapes b) et c) pour toutes les positions des recepteurs pour 
20 lesquelles on veut un resultat de demigration ; et 

e) on repete les etapes a) a d) pour toutes les valeurs prises par le parametre 
(px) pour lesquelles on a pointe un evenement sur les donnees migrees correspondantes. 

2) Methode pour determiner les temps de trajet en reflexion d'evenements 
sismiques pointes sur des enregistrements en 3D de donnees sismiques correspondant a des 

25 signaux captes par des recepteurs sismiques repartis suivant plusieurs lignes d'acquisition 
en reponse a remission dans le sous-sol d'ondes depuis un ou plusieurs points-source, ces 
donnees ayant ete au prealable transformees en une superposition de donnees en ondes 
cylindriques migrees et interpretees, caracterisee en ce que : 
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a) on definit un vecteur lenteur ( p ) dont la composante (px) suivant une 
direction parallele aux lignes d' acquisition, definit la pente des ondes cylindriques 
associees aux differentes lignes d' acquisition ; 

b) pour un recepteur sismique donnee situe a une abscisse ( .r^ ) sur une ligne 

5 d' acquisition, on recherche I'abscisse ) d'un point source sur la ligne d'acquisition telle 

qu'un rayon partant du recepteur sismique et se reflechissant sur un evenement pointe, 

— > 

emerge au dit point source, avec un vecteur lenteur ( p ) dont la composante suivant la 
ligne d'acquisition est (-px) ; 

c) on determine un temps de trajet ( t' (x^ ) ) en ajoutant a la valeur du temps de 
10 trajet le long du dit rayon, un temps egal au produit de la pente de Tonde cylindrique par 

I'abscisse du point source ; 

d) on repete les etapes b) et c) pour toutes les positions des recepteurs pour 
lesquelles on veut un resultat de demigration ; et 

e) on repete les etapes a) a d) pour toutes les lignes d'acquisition pour 
15 lesquelles on veut un resultat de demigration ; 

f) on repete les etapes a) a e) pour toutes les valeurs prises par le parametre 
(Px) pour lesquelles on a pointe un evenement dans les donnees migrees correspondantes, 

3) Methode selon Tune des revendications 1 ou 2, caracterisee en ce que Ton utilise 
les temps de trajet associes aux evenements en ondes cylindriques, pour mettre en ceuvre 

20 une technique de tomographic de reflexion par ondes cylindriques et calculer une 
distribution des vitesses dans le milieu. 

4) Methode selon Tune des revendications 1 ou 2, caracterisee en ce que Ton 
convertit les temps de trajet associes aux evenements en ondes cylindriques, en temps de 
trajet associes aux evenements par points de tir, en exploitant les correspondances entre les 

25 donnees en ondes cylindriques et les donnees en points de tir. 

5) Methode selon la revendication 4, caracterisee en ce que, les donnees sismiques 
etant organisees par lignes d'acquisition, on transforme les temps de trajet par tir synthetise 
en des temps de trajet associes aux couples source - recepteurs , cette transformation se 




faisant par extrapolation suivant la direction des lignes d'acquisition, du temps de trajet 
synthetise pour le couple source - recepteur le plus proche du couple source - recepteur du 
dispositifd' acquisition, que Ton veut documenter en temps de trajet. 

6) Methode selon la revendication 4, caracterisee en ce que, les donnees sismiques 
5 etant a azimut constant, on transforme les temps de trajet par tir synthetise, en des temps de 

trajet associes aux couples source - recepteurs , cette transformation se faisant par une 
double extrapolation, la premiere suivant la direction definie par T azimut, la deuxieme 
suivant la direction orthogonale a la direction precedente. 

7) Application de la methode selon Tune des revendications 4 a 6, a la 
10 determination du modele de vitesse, en appliquant par exemple une tomographie de 

reflexion aux temps de trajet obtenu. 
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